Decoherence in an infinite range Heisenberg model 



A. Dey*, M. Q. Lone* and S. Yarlagadda 

TCMP Div., 1/AF Salt Lake, Saha Institute of Nuclear physics, Kolkata, India. 
E-mail: y.sudhakar@saha.ac. in 

PACS numbers: 03.65.Yz, 03.67.Pp, 75.10.Jm, 05.70.Ln 

Abstract. We study decoherence in an infinite range Heisenberg model (IRHM) 
in the two situations where the system is coupled to a bath of either local optical 
phonons or global optical phonons. Using a non-perturbative treatment, we derive 
an effective Hamiltonian that is valid in the regime of strong spin-phonon coupling 
under non-adiabatic conditions. It is shown that the effective Hamiltonian commutes 
with the IRHM and thus has the same eigenstates as the IRHM. By analyzing the 
dynamics of the system using a quantum master equation approach, we show that 
the quantum states of the IRHM system do not decohere under Markovian dynamics 
when the spins interact with local phonons. For interactions with global phonons, 
the off-diagonal matrix elements of the system's reduced density matrix, obtained for 
non-Markovian dynamics, do not indicate decoherence only when states with the same 
Si^ (i.e., eigenvalue for the z-component of the total spin) are considered. 
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1. Introduction 

Quantum information processing heavily relies on a precious and fragile resource, 
namely, quantum entanglement pQ. The fragility of entanglement is due to the 
coupling between a quantum system and its environment; such a coupling leads 
to decoherence, the process by which information is degraded. Decoherence is the 
fundamental mechanism by which fragile superpositions are destroyed thereby producing 
a quantum to classical transition [2J |3]. In fact, decoherence is one of the main 
obstacles for the preparation, observation, and implementation of multi-qubit entangled 
states. The intensive work on quantum information and computing in recent years 
has tremendously increased the interest in exploring and controlling decoherence effects 

ilia ei Ei ei is nam. 

Since coupling of a quantum system to the environment and the concomitant 
entanglement fragility are ubiquitous [U 12], it is imperative that progress be made 
in minimizing decoherence. Decoherence free states prevent the loss of information due 
to destructive environmental interactions and thus circumvent the need for stabilization 
methods for quantum computation and quantum information. In the past decoherence- 
free-subspace (DFS) (121 [13] has been shown to exist in the Hilbert space of a model 
where all qubits of the quantum system are coupled to a common environment with 
equal strength. A DFS is a subspace which is invariant under the action of the system 
Hamiltonian; furthermore, the subspace is spanned by degenerate eigenvectors of the 
system operators coupling to the environment [TH [15]. Alternately, decoherence can 
also be suppressed through quantum control strategies [T6l fT71 ITS] . 

Although the theory of decoherence has undergone major advances [21 [3], yet, there 
exist many definitions of decoherence [19] . For the analysis in this paper, we choose the 
most commonly used definition of decoherence: Loss of off-diagonal elements in the 
system's reduced density matrix. In general, a many-qubit (i.e., many-spin) system can 
have distance dependent interaction. The two limiting cases for interaction are spin 
interactions that are independent of distance and spin chains with nearest-neighbor 
interactions only. In this work we consider the extreme case of distance independent 
interaction among the spins, i.e., the IRHM. The objective of this paper is to study 
the decoherence phenomenon, due to coupling of spins of the IRHM to optical phonons, 
in the two extreme cases of the spins being independently coupled to different baths; 
and all the spins being collectively coupled to the same environment. We employ 
the analytically simpler frame of reference of hard-core-bosons (HCBs) rather than that 
of spins so that the single particle excitation spectrum can be easily obtained and 
exploited; we show that the effective Hamiltonian even in higher order (i.e., greater 
than second order) perturbation theory retains the same eigenstates as the IRHM when 
the spins are coupled to local phonons. Furthermore, decoherence is studied using the 
quantum master equation approach [2D]- Our dynamical analysis shows that the system 
coupled to local phonons does not decohere when Markov processes are considered; 
whereas for global phonons, even for non-Markovian dynamics, there is no decoherence 
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when eigenstates with the same eigenvalue S T (i.e., z-component of the total spin) are 
considered. 

The rest of the paper is organized as follows: In section 2, we introduce the 
IRHM Hamiltonian and describe its eigenstates and eigenenergies. In section 3, we 
study decoherence under strong coupling with local optical phonons, and show that the 
effective Hamiltonian thus obtained retains the same eigenstates as Hirhm- in section 4, 
we use the master equation approach and show that the system does not decohere under 
local and global couplings. Next, in section 5, we give our conclusions and make some 
general remarks regarding the wider context of our results. The paper also contains 
an appendix where we derive the third order perturbation contribution to our effective 
Hamiltonian (H e ff) and show that the eigenstates of the IRHM Hamiltonian are retained 
by our H eff . 



2. Infinite Range Heisenberg Model 



We begin by introducing the IRHM whose decoherence will be studied when the system 
is coupled to either local or global optical phonons. The IRHM is defined as: 



Hirum = J ~/ \ Si.Sj + (A — l)S?Sj 



i,j>i 




(1) 



where J > 0, A > 0, and we are considering only S 



commutes with both S Total 



1/2 spins. We note that Hibmm 

J2i s ?) and (Ei<Si) (= S T otai)- In equation CQ), it is 
understood that J = J*/(N — 1) (where J* is a finite quantity) so that the energy per 
site remains finite as iV — >■ oo. The eigenstates of -^irhm are characterized by St (i.e., 
the total spin eigenvalue) and S T (or the eigenvalue of the z-component of the total spin 
; the eigenenergies of these eigenstates are 
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(2) 



The ground state corresponds to S T = and St = which is rotationally invariant. 



3. Effective Hamiltonian for IRHM spins coupled to local optical phonons 

The real quantum computer will not be free from noise and thus the entangled states 
have a tendency to undergo decoherence. To study decoherence due to phonons, 
we consider interaction with optical phonons such as would be encountered when 
considering transition metal oxides. We will now derive an effective Hamiltonian, 
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when the spins of the IRHM are coupled to local optical phonons, and show that the 
eigenstates of the effective Hamiltonian are the same as those of -£/irhm- 
The total Hamiltonian is given by 

H T = #irhm + gu ^2 S i( a l + a i) + u ^2 a ^ ah ^ 



where a is the phonon destruction operator [21], uj is the optical phonon frequency, 
and g is the coupling strength. Now, we make the connection that the spin operators 
can be expressed in terms of HCB creation and destruction operators and b, i.e., 
fet = S + , b = S~, and tfb = S z + 0.5. We then observe that conservation of S^ otal 
implies conservation of total number of HCB. The total Hamiltonian is then given by 



h = jj2i(°M b j + H - c -) + A o* - °- 5 )K' - °- 5 )] 

i,j>i 



J 3 

where n.j = bjbj. Subsequently, we perform the well-known Lang-Firsov (LF) 
transformation [221 123] on this Hamiltonian. Under the LF transformation given by 
e s He~ s = H + Hj with S = —g — |)( a « — a \)y ^ ne operators bj and dj transform 

like fermions and bosons; this is due to the interesting commutation properties of HCB 
given below: 

M,-] =[6i,6t] = 0, fori^i, 

{6 i? 6l} = l. (5) 
Next, the unperturbed Hamiltonian Hq is expressed as [23] 

Hq = H s + H env , (6) 
where we identify H s as the system Hamiltonian 
H s = J [(0.5e- ff2 6j^ + H.c.) 

i,j>i 

+ A(n i -0.5)(n i -0.5)] J (7) 



and H env as the Hamiltonian of the environment 

H env = ^y^gtqj. (8) 



On the other hand, the interaction Hi which we will treat as perturbation is given by 
Hi = jJ2[ Q - 5e ~ 92 b\bj]{Sf S 1 ! - 1} + H.c, (9) 

i,j>i 

where S± = exp[±p(aj— a^-)]. In the transformed frame, the system Hamiltonian depicts 
that all the HCBs are coupled to the same phononic mean-field. Thus, the unperturbed 
Hamiltonian H comprises of the system Hamiltonian H s representing HCBs with the 
same reduced hopping term 0.5Je _s and the environment Hamiltonian H env involving 
displaced bath oscillators corresponding to local distortions. Here it should be pointed 



Decoherence in an infinite range Heisenberg model 



5 



out that both the interaction of the HCB with the mean-field as well as the local 
polaronic distortions in the bath oscillators involve controlled degrees of freedom. Now, 
the system Hamiltonian H s can be expressed as 

H s = #irhm + (H s — Hirhm) (10) 

When we change the Hamiltonian from Hjrbm to H s by adiabatically turning on the 
perturbation (H s — -£/trhm), the resulting state of the system is still obtainable from 
that of i^iRHM by using unitary Hamiltonian dynamics and is thus predictable based 
on a knowledge of the coupling parameter g [21] . Thus no irreversibility is involved in 
going from -Htrhm to H s . On the other hand, perturbation Hi pertains to the interaction 
of HCBs with local deviations from the phononic mean-field; the interaction term Hj 
represents numerous or uncontrolled environmental degrees of freedom and thus has the 
potential for producing decoherence. Furthermore, it is of interest to note that the 
interaction term is weak in the transformed frame compared to the interaction in the 
original frame; thus one can perform perturbation theory with the interaction term. 

We represent the eigenstates of the unperturbed Hamiltonian H as \n,m) = 
\n) s ®\m) p h with the corresponding eigenenergies E n>m = E^ + E^ 1 ; \n) s is the eigenstate 
of the system with eigenenergy E s n while \m) p h is the eigenstate for the environment 
with eigenenergy Ef£. Henceforth, for brevity, we will use u m = E^. On observing 
that (0, 0|i7j|0, 0) = ( i.e., the ground state expectation value of the deviations is 
zero), we obtain the next relevant second-order perturbation term 

(0,0|#j|n,m) (n,m\Hi\0,0) 
Eon — E nm 



f(2) = y (U,U|^j|n,m)(ra,m|Jjj|U,U) 



For strong coupling (g > 1) and non-adiabatic (J*/ui < 1) conditions, on noting that 
uj m — ujq = u) m is a positive integral multiple of oo and that E s n — E^ ~ J*e~ 9 uj 
(as shown in the next section), we get the following second-order term [25] using 
Schrieffer- Wolff (SW) transformation (as elaborated in Appendix A of references 
and 



ph (0\Hj\m) ph ph (m|iJ 7 |0) ph 



m 



E 

i,j>i 



where 



[(0.5J| 2) ^ + H.c.) 

-0.5j{ 2) {ni(l - n 3 ) + 72,(1 - m)}] , (12) 
Jf s -(Ar-2)/ I ( 9 )^~-( W -2)^, (13) 

/.f'-PAW + Ato)]^-^, (14) 

with Mg) = Er.i9 2 "/(« ! «) and Ml) = E"=iES.iS 2( " + "' , /[n!m!(n + m)]. The 
effective Hamiltonian H s + is a low energy Hamiltonian obtained by the canonical 
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SW transformation decoupling the low-energy and the high-energy subspaces; 

2 

this decoupling is a consequence of J*e~ 9 <C u. We now make the important 
observation that the effective Hamiltonian H s + H^ 2 \ when expressed in terms of spins, 
has the following form: 

[MSfS? + S»SJ) + J,„ g SfS«] , (15) 



E 

i,j>i 



and thus has eigenstates identical to those of the original Hamiltonian -£/irhm hi equation 
([1]) because £\ j^S- : Sj) and -£/irhm commute. On carrying out higher order (i.e., 
beyond second order) perturbation theory (as discussed in Appendix A), and expressing 
the results in the spin language, we still get an effective Hamiltonian H e ff of the following 
form that has the same eigenstates as the IRHM. 



j 

i,j>i . 



i k 

where J xy and J z are functions of the S^ otal (= J2k $k ) operator. The small parameter 
of our perturbation theory, for a small N system, is J/(gu) [see reference [30J for 
details]; whereas for a large N, the small parameter is J*/ (g 2 oo) [see reference [31] for an 
explanation] . It is the infinite range of the Heisenberg model that enables the eigenstates 
of the system to remain unchanged. Next, we study decoherence in a dynamical context 
and gain more insight into how the states of our -^irhm can be decoherence free. 



4. Dynamical evolution of the system 

In this section, we will study decoherence in the system from the dynamical 
perspective. We will discuss the dynamics of an open quantum system, described by 
the Hamiltonian IRHM, using master equation approach. Our quantum system is open 
because it is coupled to another quantum system, i.e., a bath or environment [34]. In 
our case, IRHM is coupled to a bath of either local optical phonons [see equation ([3])] or 
global optical phonons [see equation fl38l) ]. As a consequence of the system-environment 
coupling, the state of the system may change. This interaction may lead to certain 
system-environment correlations such that the resulting state of the system may no 
longer be represented in terms of unitary Hamiltonian dynamics. The dynamics of the 
system, described by the reduced density matrix p s (t) at time t, is obtained from the 
density matrix pr{t) of the total system by taking the partial trace over the degrees of 
freedom of the environment: 

p s (t) = Tr R [p T (t)} = Tr R [U(t)p T (0)U\t)} , (17) 

where U (t) represents the time-evolution operator of the total system. Now it is evident 
from the above equation that we need first to determine the dynamics of the total system 
which is a difficult task in most of the cases. By contrast, master equation approach 
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conveniently and directly yields the time evolution of the reduced density matrix of the 
system interacting with an environment. This approach relieves us from the need of 
having to first determine the dynamics of the total system-environment combination 
and then to trace out the degrees of freedom of the environment. 

4-1. Decoherence due to Local Optical Phonons: 

We begin this sub-section by considering the following Hamiltonian: 

H = H + H h (18) 

where Hq is the system-environment Hamiltonian given by equation ([6]) and Hj 
represents the interaction Hamiltonian given by equation Q). It is convenient and 
simple to derive the quantum master equation in the interaction picture. Thus our 
starting point is the interaction picture von Neumann equation for the total density 
operator prif) 

*^ = -i[H I (t),jh{t)\, (19) 
where Hi(t) = ^ H ot^ ie -iU t anc j p T (t) = e lHot p T {t)e~ %Hot are the interaction 
Hamiltonian and the total system density matrix operators (respectively) expressed 
in the interaction picture. Re-expressing the above equation in integral form yields 

~p T (t) = p T (0) - 1 [ drlH^rPTir)]. (20) 



Jo 

Nowadays there is considerable interest in systems with initial correlation with 
the environment [321 [33]; however, for simplicity, let us suppose that the initial 
state of the total system is a factorized state given as />r(0) = p s (0) (8> Ro with 
Ro = Yin \ n )ph P h(n\e~P Un / Z being the initial thermal density matrix operator of the 
environment and (3 = -j^f] furthermore, Z = Y in e~ l3uJn defines the partition function 
of the environment. With this assumption, we substitute equation (120]) inside the 
commutator of equation (Tl9l) and then take the trace over the environmental degrees of 
freedom to obtain the following equation: 

< ^ L = -iTr R [H I (t) t p.(0)®R o ] 

^ drTrniHj^^Hj^^Tir)}}. (21) 







The above equation still contains the total density matrix Pt{t); In order to evaluate 
it, we rely on an approximation known as the Born approximation. This approximation 
assumes that the environment degrees of freedom are large and thus the effect on 
the environment due to the system is negligibly small for a weak system-environment 
coupling. As a consequence, we write Pt(t) = Ps{t)®Ro+0(Hi) within the second order 
perturbation in system-environment interaction pH ESI [36l [371 EHl EH1 HO, HU H21 133] . 
Therefore we can write the equation (|2"T|) in time-local form as 

-zTr R [Hj(t), Ps (0)^R o ] 
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- [ t dTTr R [H I (t),[H I (T),p,(t)®R ]l (22) 
Jo 

We note here that, for obtaining the non-Markovian time-convolutionless master 
equation ff22l) . we replaced p s (r) with p s (t). This replacement is equivalent to obtaining 
a time-convolutionless master equation perturbatively up to only second order in the 
interaction Hamiltonian using the time-convolutionless projection operator technique 
[34"t |39| HQ]. It has been shown in a number of cases that time- local approach works 
better than time-nonlocal approach [341 [371 HU H2, HI] . Now we will consider the second 
order time-convolutionless master equation (122]) with the time variable r replaced by 
(t-r). 

dp s (t) 



dt 

ft 



[ dTTr R [H I (t),[H I (t-T),p s (t)®R }}. (23) 
Jo 

Next, we will study the Markovian dynamics of the system. To this end we assume 
that the correlation time scale r c for the environmental fluctuations is negligibly small 
compared to the relaxation time scale t s for the system, i.e., r c <C t s . This time scale 
assumption is motivated by the condition J*e~ g <C uj already mentioned in section 3. 
The Markov approximation (r c <C t s ) allows us to set the upper limit of the integral to 
oo in equation f l23|) . Thus we obtain the second order time-convolutionless Markovian 
master equation ( l24l) : 

df ^= -iTr R [Hr(t),p s (0)®R o ] 

POO 

- / drTr R [H I (t),[H I {t-r),p s (t)®R ]\. (24) 
Jo 

Defining {|?i)p/i} as the basis set for phonons, therefore, we can write the master equation 
as: 

P h(n\[Hi(t),p s (0)®R o ]\n) ph 



poo r 

dr ph (n\Hj(t)Hj(t - r)p s (t) <g> R \n) ph 
... Jo l 



- ^(nlHj^Psit) g> RoH^t - r)\n) ph 

- - r)p s (t) g> RoHj^lnjph 

+ P h(n\p s (t) ® RoH T (t - r)^j(t)|n)ph] . (25) 

In order to simplify the above master equation, we need to evaluate the time 
evolution of the operators involved in Hi. Considering the second term in the equation 
(T25D, yields 

pft (n|F/(f).F/j(£ - r)p s (t) <g> R \n) ph 
= Y, elHst pMHi\m) ph e- iH ^ e ^^^H/T 7 |n> 1 *e- iff -<*- T )^(t)-^- e < ^- w »J T .(26) 

m 
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J2k eikrj *>t and 



In momentum space, we express HCB operators as. ^ — ^ k o - u k 
bj = -7= e~ lkrj bk] then, it is important to note that the hopping term in the system 
Hamiltonian can be written as: 



0.5jJ2( e ~ g2b l b J + H.c.) = 0.5 Je 



i,j>i 



N , 



e" 5 n - 0.5Je~ 9 N p 



(27) 



where we used J = J*/(N — 1), N p = ^2 k b\.bk and no = b\b a (i.e., the particle number 
in momentum k = state). Here it should be mentioned that using HCBs instead of 
spins has enabled us to obtain (with ease) the excitation spectrum which is crucial 
for the analysis given below. Let {|<?) s } denote the complete set of energy eigenstates 
(with eigenenergies of the system Hamiltonian H s ; then we can write: 



AH s t 



H ie 



-iH s t 



l,j>l q,q> 



N 2^ °k°V e 

+ H.c, 



■* B ' t \4).Jtf\{s$s!> 



which implies 

9,9' 

(29) 

where \E* — Et\ = 0-5J*(j^j)e~ g2 or . Here we have taken the total number of 
HCBs to be conserved; then, only the hopping term in H s will contribute to the particle 
excitation energy. Substituting equation ( J29|) in equation ( 1261) . we get 

pMH^H^t - r)p s {t) ® R \n) ph 

m q,q',q" 

x p s (t)^— e i(w ' l " Wm)T . 

(30) 

Thus under the assumption of J*e~ 9 << uj, it follows that \u n — oo m \ » \E* — Et\ 
and \u n — uj m \ » \Et — Et,\; hence in equation (1301) . we can take e q i = 1 
and e l ^ E q'~ E q" S) ( t ~ T ^ = i which implies that we do not get terms producing decay. The 
resultant equation is 



h (n\Hi(t)H!(t - r)p s (t) ® R \n) ph = ^ ph {n\H ^m) ph ph (m\H >|n) pft p B (t) 



e 



-e 



i(a) n — a; rrl )T 



(31) 
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Carrying out the same analysis on the remaining (i.e., third, fourth, and fifth) terms 
in the master equation, we write equation (1251) as: 



dt 



P h(n\[Hi(t),p a (0)®Ro]\n) 



ph 



E/ * 



ph 



[n\Hj\m} ph ph^H^njph p s (t) e 



i(u n — w m )r 



ph 



ph 



[n\H!\m) ph p s (t) ph (m\H I \n) i 



ph " 



[n\Hi\m)ph, p s {t) v h(m\Hi\n)ph- 



+ Ps(t) ph^H^mjph ph {m\Hi\n) 



-fan 



D i(ul n -U) m )T 



-i(u„-U) m )T 



ph ' 



-i{w n —U) m )T 



• (32) 



Next, we evaluate the first term in the above equation and show that it is zero at T = 0. 
We observe that 



Tr R [Hj(t)R } = P h(n\Hi(t)R \n) 



ph 

= 0.5Je^ 2 [e iHst bjb je - iH ^ ph {0\{S l f S l l - l}\0) ph 

= 0. (33) 

Thus, we have p h{n\ [Hj(t), p s (0) <S> R ]\n) p h = and the master equation at T = 
simplifies as: 

^ir= ~E /" ' dT[\ ph (0\Hj\m) ph \ 2 p„{t)e**" + p s (t) |^<0|ff / |m) pfc |V u » T ] 

m 

poo 

+ J2 dT [phHHmphPsitUiOlHjl^p^ 

n Jo 

+ ph (n\Hj\ 0) ph p s (t) ph (0 1 if j | n) ph e- WnT ] 
dr e^^ T | p/l (0|/f 7 |n) p ,| 2 p s (t) 

POO 

+ / dr e l ^ + ^ T p s (t) IphiOlHj^phl 2 
Jo 



E 



dr e""- T ph (n\H,\0) rh p,(t) ph {0\H,\n) 



ph 



(34) 



Now, we know J™ dre luJ " T oc S(u n ). Therefore, on using this relation and the fact that 
p h(0|iff|0)ph = 0, the third term in equation ( 13~4"|) vanishes; hence, we get 



dp s {t) 
dt 



E 



The term Y^ r 



\ P h{0\Hi\n) ph \ 



\ph 



(0|if/|n) 



ph | 



OJr, 



■p a {t) -p a (t) 



I ph 



(0\Hi\n) 



pii I 



1 2 n 



(35) 



corresponds to the effective Hamiltonian if ( 2 ) in second 



order perturbation and commutes with ifg (see section 3). Let \n} s be the simultaneous 
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eigenstate for and H s with eigenvalues E n and E n , respectively. Then, from the 
above equation we get: 

s (n\p s (t)\m) s = s ( n \p s (0)\m) s , (36) 

which implies that 

,(n\Ps(t)\m). = e -^- E ^ s (n\p s (0)\m) s , (37) 

(2) 

where E n = E n + E n . Thus we see from the above equation that there is only a 
phase shift but no decoherence! Since the matrix elements of an operator are invariant 
under canonical transformation, it should be clear that no loss in off-diagonal density 
matrix elements (i.e., no decoherence) in the LF transformed frame of reference implies 
no loss in off-diagonal density matrix elements (i.e., no decoherence) in the original 
untransformed frame of reference. Although the HCB's in the original frame of reference 
form polarons and are thus entangled with the environment, nevertheless no decoherence 
results. For greater clarity, the form of s ( n \ps{t)\ m ) s in the original frame of reference 
and its associated non-decoherence is discussed in Appendix B for a special two-spin 

2 

case of IRHM. Thus, up to second order in perturbation, the assumption J*e~ 9 << oj, 
the infinite range of the Heisenberg model, and the Markov approximation (r c <C r s ) 
together have ensured that the system, with a fixed Sj>, does not decohere. 

While the above analysis is valid in the regime ksT /u « 1, the finite temperature 
case ksT/uj > 1 needs additional extensive considerations and will be dealt with 
elsewhere [15]. 



4-2. Decoherence due to Global Phonons: 

We will now analyze decoherence due to interaction of the spin system with global 
phonons. To this end, we consider the following total Hamiltonian where all qubits of 
our IRHM interact identically with the environment. 

if Tot = H s + ^ 5? ^ w *(^fc4 + 9k a k) 

i k 

+ ) j u) k ala k , (38) 
k 

where H s = Hirhm is the Hamiltonian of the system. (Here, for global phonons, since we 
do not use LF transformation, we define H s as the untransformed system Hamiltonian.) 
Since the z-component of the total spin S^ otal (and thus the interaction Hamiltonian) 
commutes with if irhm > the eigenstates having same eigenvalue S* T constitute a DFS. To 
study the case when S* T is not conserved and to obtain the form of the reduced density 
matrix p s {t), we study the dynamics of the system through the following non-Markovian 
master equation [35] : 

dp s (t) 



dt 



-i[H s ,p s (t)] + F(t)[Lp s (t),L] 
+ F*(t)[L,p s (t)L], (39) 
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where L is the system operator that couples with the bath and satisfies the constraint 
[L,H S ] = 0. For the total Hamiltonian in equation fl38j) . L = ^i^t = S^ otal . Also, 
F(t) = Jq a(t — s)ds where a(t — s)— r)(t — s) + iv{t — s) is the bath correlation function 
at temperature T with 




v{t-s)= -J2 \9k?^l sin[w fc (i - s)}. (40) 

k 

The function F(t) governs the non-Markovian dynamical features of the system. 

Let {|n) s } be the eigen basis in which both the operators S^ otal and H s are 
simultaneously diagonalized. Upon solving the master equation explicitly we get [4"6] : 

s (n\p s (t)\m) s 

= exp (-z l(K - E s Jt + {(S^f ~ (SU 2 } Y(t)]) 
x exp [- {S z Tn - S z Tm fX{t)\ s {n\p s (0)\m)s, (41) 

where and S^ n are defined through H s \n) s = E^\n) s and ^ 4 Sf\n) s = S^ n \n) s . 
Furthermore, X(t) = f* F R (s)ds and Y(t) = / Q * F/(s)ds with F R {t)+iF T {t) = F{t). This 
implies that, when states \m) s and \n) s have the same z-component of the total Spill Srp 
( i.e., S^ n = S^ m ), the matrix elements s (n\p s (t)\m) s = s (n\p s (0)\m) s exp[-i(E^ -E s m )t] 
display a decoherence free behaviour. But when is not conserved, the off-diagonal 
matrix elements will diminish in general, i.e., the system undergoes decoherence. In the 
language of HCBs, the eigenstates of the system with a fixed number of HCBs makeup 
a DFS. Furthermore, the entanglement entropy of the system will remain unaltered 
since the density matrix evolves unitarily. In future, using the above framework, we 
will consider the interesting case of dynamical evolution and decoherence of states with 
different values. 

5. Discussion and Conclusions 

In conclusion, we have shown that the eigenstates of H e ff are the same as those of 
-f^iRHM an d for Markov processes they are decoherence free under the coupling of the 
system to local optical phonons. For global optical phonons (i.e., when all the qubits 
are exposed to the same collective noise) the eigenstates with the same form a DFS. 
A DFS is expected in the global phonon case because the Hamiltonian of the system 
commutes with Sf. But the important point is that, even in the local phonon case, 
it is still possible to fully preserve coherence for the composite particle (i.e., polaronic 
HCB) system with a fixed number of particles. More specifically, for local phonons, 
s( n \ps{t)\ m ) s differs from s {n\p s (0)\m) s only by a phase factor and s (n|p s (0)|m) s can be 
obtained from s (ra|piRHM| m )s (density matrix element of IRHM) by an exact unitary 
evolution [21]. Later, we will analyze the non-Markov processes and see how the 
resultant dynamics deviates from the Markovian dynamics. 
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Earlier, a new type of resonating valence bond (RVB) states |47] were constructed 
for four and six spins using homogenized linear superposition of the St = states of 
-f^iRHMi these RVB states have a high bipartite entanglement. The decoherence analysis 
in this paper is also applicable to these new RVB states which are groundstates of our 
#irhm- Our RVB states are constructed using valence bond (VB) states which are 
St = states. VB states are built from singlet states between pairs of spins. A general 
VB state is defined as: 

|*)vb = \® iujl ) <8> \®i 2 , j2 ) <8> ... (g> \$ iM j M ), 

where \&i k> j k ) = ^ (If)iJ ~ \)j k ~ I ~ denotes the singlet dimer connecting 

a site z'fc in sub-lattice A with a site j k in sub-lattice B. Examples of our RVB states 
(that are constructed from spins 1, 2, 3, 4,... arranged sequentially on the vertices of a 
regular polygon and that have high bipartite entanglement) are l^)^ given below for 
four spins: 

|^ 4 ) rvb = cj 3 (|*i, 2 ) ® |*3,4)) + wf (|$i j4 ) ® |*2,3)), (42) 
where co 3 (= e i2 ^/ 3 ) is a cube root of unity; and |*6)rvb given below for six spins: 

|* 6 )rvb = ^4(1*1,2) ® |*3, 6 ) ® 1*4,5)) 
+ W 4 2 (|$2, 3 ) ® |*1,4) ® |*5, 6 )) 

+ ^(|*1, 6 ) ® 1*2,5) ® 1*3,4)), (43) 

where (= e l27r//4 ) is a fourth root of unity. 

Before closing we will make a few general remarks. Firstly, the -£^irhm model of 
this paper, deals with the extreme case of distance independent interaction among the 
spins. On the other extreme end, if one were to consider a nearest-neighbor interaction 
anisotropic Heisenberg chain [of the type ^2i{J±(SfSf +1 + SfSf +1 ) + J\\S?Sf +l } where 
J_i_ > J\\ > 0] with a strong coupling to local phonons [introducing the additional terms 
gu Sf(a\ + a,i) + oj ^ a\ai\, then the spin system undergoes a Luttinger liquid to a 
spin-density-wave transition upon turning on the spin-phonon interaction and decoheres 
[4*8] . In general, distance-dependent-interaction in spin Hamiltonians will fall somewhere 
in between the above two extreme cases. 

Next, our decoherence analysis for local optical phonons will continue to be valid 
even for the more general optical phonon terms given below: 

^72" Yl S i\- U k(9kali + g* k a Ki )} + ^ w *4,i a M- ( 44 ) 

i,k k,i 

We also must mention that our approach cannot accommodate the acoustic phonon case 
as here the condition J*e~ 9 << Uk cannot be satisfied in the long wavelength limit. 
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Figure Al. Open loop hopping processes contributing to effective hopping term 
T% in third-order perturbation theory. Here empty circles correspond to sites with no 
particles while filled circles correspond to sites with hard-core-bosons. The numbers 
1,2, and 3 indicate the order of hopping. 




2 

(c) 



Figure A2. Closed-loop hopping processes contributing to effective interaction term 
in third-order perturbation theory Here filled (empty) circles correspond to sites 
with (without) hard-core-bosons. The numbers 1, 2, and 3 represent hopping sequence. 

Appendix A. 

In this appendix, we will show that the third-order perturbation theory also produces 
a term that has the same eigenstates as IRHM. To this end, we obtain the following 
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(c) 

Figure A3. Hopping processes (involving closed loops) contributing to effective 
hopping term Tq h in third-order perturbation theory. Filled (empty) circles represent 
occupied (unoccupied) sites. 



third-order perturbation term in the effective Hamiltonian: 

ph^H^rnjph p^mlH^ph ph {n\H >|0> ph 



AE p m h AE p n h 



(A.r 



Here AE% = uj m — ujQ. Evaluation of leads to various hopping terms and interaction 
terms. 



t n T n + t cn T ( 



1 1 

Cn 



n=l 



n=l 



EE* 



i n=l 



(A.2) 

J 3 /(g 2 u) 2 (as will be 



where t n ~ (J 3 e 9 )/(g 2 u) 2 , t cn ~ J 3 e 9 /(goj) 2 , and t> n ~ 
explained later). We will demonstrate below that is of the following form 



T(^n Jfe )6f6 i + H.c. 



E^E"*) 



(A.3) 



where T and V are functions of the total number operator Y2k n k- Since the IRHM 
commutes with the total number operator, has the same eigenstates as IRHM! 

There are six open-loop hopping processes T% depicted in figure IA11 We analyze 
them sequentially below. 

T i = E E 6 M¥»fr 
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k^i,l,j 



E^-^A)- 1 



E^- 6 ! 6 *)- 1 



j 



(iV-2)-E% 
(JV-1)-X>& 

EC 1 - 6 ! 6 *)-! 

(JV-2)-E)6l6* 



6ft 



The second hopping process in figure IA1I (b) is given by 

T * = E E 

= E (! - 6t A) E 6 M 6 « 



Ec 1 - 6 ^ 



fc^i.Z 



(^v-i)-E 6 fe 



6f6*. 
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(A.4) 



(A.5) 



The hopping process in figure lATl fc) is expressed as = J2k^i i j i 6 l 6 fc 6 ] 6 » 6 l 6 .? = 
Tl 1 . The fourth hopping process T l £ in figure [ATI (d) is obtained as follows. 

T i l = E E 6 fe 6 M 6 * 

jy^i,l,k kj^i,l 

= E (! - h %) E b l b ^ 

jj^i,l,k kj^i,l 

= T l 2 \ (A.6) 

The hopping process in figure [AT] (e) yields = Y,j 1 Li,i,kY J k^i,ib\bjb\b k b]b i = Tj\ 
We analyze below the last hopping process Tq in figure lA~TI (f). 

T 6 l = E E 6 M 6 ; 6 I 6 * 
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E E 



E 6t A - 1 



k^i,l 



E ^ - 1 



E 6 M*« 

. 3 

E^A-2 



(A.7) 



We will now deal with closed-loop hopping processes such as those in figure IA21 
These lead to effective interactions. The process V[ in figure IA2I (a), obtained from 
figure [ATI (a) by setting I = i, is given as follows. 

vt = EE & JMfe^ 

= E( 1 -^)E( 1 -%) 6 ^ 



E(! - 6 &) - 1 



w-E%- 



(^-i)-E 6 ^ 



(A.8) 



Next, the hopping process corresponding to closed loop in figure lA"2l (b) is obtained 
from figure IA1I (c) by taking / = %. 



D 1 - ^ 



6l6i. 



(A.9) 



Lastly, the hopping V3 [depicted by the closed loop in figure IA2I (c)] is obtained from 
figure IA1I (e) by setting I = i. 



v * = EE 6 & 6 Mfc 



j=£i,k k=£i 

E(i-^)E^A 



k^i 
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(A.10) 

Finally, we consider figures IA~3l (a) . (b), and (c) which deal with effective hopping 
terms T^ n involving closed loops. The effective hopping term T^, corresponding to 
figure lA3l (a), is obtained by setting k = i in figure lATl (a): 



■ci 



E 



(iV-2)-E% 



(iV-l)-^^ 



b\h 



b\b t . 



To obtain the effective hopping term T^ 2 corresponding to figure 
in figure IA1I (a): 



(A.ll) 

(b), we take j = I 



T, 



C2 



E h W h Mh 

£(1 

k^i,l 



(N — 2) 



(N 



k^l 

£fe 



b\h 



b)h 



rpli 

J ci- 



The effective hopping term Tq 3 depicted in figure IA3I (c) [upon setting k 
in figure IA1I (a)] is given by 



Tea = b\b % b\hb\b t = b%. 



(A.12) 
i and j = I 

(A.13) 

Thus we have shown that contains effective hopping terms £\ i >i [T(^2 k n^bjbi 
+H.c] and effective interaction terms V^Ylik n k) n i)- Since T and V are functions of 
the total number operator, and IRHM have the same eigenstates. These arguments 
can be extended to even higher-order perturbation theory to show that the effective 
Hamiltonian (after taking all orders of perturbation into account) will give the same 
eigenstates as IRHM! 

We will now explain the expressions for the coefficients t n , v n , and t cn in equation 
(IA.2I) . obtained from third-order perturbation theory, using typical schematic diagrams 
shown in figure IA4I [for details of corresponding diagrams and analysis in second order 
perturbation, see reference [26]]. We consider two distinct time scales associated with 
hopping processes between two sites: (i) ~ l/{Je~ g2 ) corresponding to either full 
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Figure A4. Schematic diagrams (a), (b), and (c), corresponding to the hopping 
processes depicted in figure IA1I (a), figure IA2I (a), and figure IA3I (a), respectively, 
yield coefficients t n , v n , and t cn , respectively. The intermediate states give the typical 
dominant contributions. Here empty circles correspond to empty sites, while filled 
circles indicate particle positions. Parabolic curve at a site depicts full distortion at 
that site with corresponding energy —g 2 Lo (+g 2 u)) if the hard-core-boson is present 
(absent) at that site. 
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distortion at a site to form a small polaronic potential well (of energy —g 2 oo) or full 
relaxation from the small polaronic distortion and (ii) ~ 1/ J related to negligible 
distortion/relaxation at a site. The coefficient t n corresponds to the typical dominant 
distortion processes shown schematically in figure IA4I (a) with the pertinent typical 
hopping processes being depicted in figure IA1I (a). In figure IA4I (a), after the HCB 
hops away from the initial site, the intermediate states have the same distortion as the 
initial state. Next, when the HCB hops to its final site there is a distortion at this 
final site with a concomitant relaxation at the initial site. Hence the contribution to 
the coefficient t n becomes J/(2g 2 u) x J/(2g 2 u) x Je~ gZ ~ J 3 e~ g2 /(g 2 u) 2 . As regards 
coefficient v n , it can be deduced based on the typical dominant hopping-cum-distortion 
processes depicted in figure IA4I (b) which typifies the hopping processes in figure IA2I 
(a). In figure IA4I (b), when the particle hops to different sites and reaches finally the 
initial site, there is no change in distortion at any site. Hence v n can be estimated 
to be J/(2g 2 u) x J/(2g 2 u) x J~ J 3 /(g 2 u) 2 . Lastly, we obtain the coefficient t cn by 
considering the typical dominant diagram in figure IA4I (c) corresponding to the typical 
process in figure EH] (a). In figure EH (c), where the first intermediate state depicts the 
particle hopping but leaving the distortion unchanged, we get a contribution J/(2g 2 u); 
for the next intermediate state, where the HCB returns to the initial site, the initial 
site has to undergo a slight relaxation (involving absorbing a phonon so as to yield a 
non-zero denominator in the perturbation theory) leading to the contribution J/ou; and 
lastly, when the HCB hops to the final site, there is a distortion at the final site with a 
simultaneous relaxation at the initial site thereby producing a contribution Je~ 9 . Thus 
we calculate t cn to be J/(2g 2 u) x J/cu x Je~ g2 ~ J 3 e~ g2 /(gu) 2 [30] . 



Appendix B. 



In equation ( 157)) . the matrix element s (ro|p s (t)|m) s can be written as 



,{n\p a (t)\m) s = s (n\ 



P h(n\p T {t)\n) ph 



\m) t 



n 



Y P h{n\e s p° T {t)e s \n) ph \m) f 



(B.l) 



where p^{t) is the total density matrix in the original frame of reference . Now, we 
illustrate this quantity by considering the simple two-spin (i.e., N=2) case of the IRHM. 
The singlet state 4g(| t4) - | 4t)) and the triplet state 4j(| 14) + I It}) are the = 
eigenstates of the two-qubit IRHM Hamiltonian; in HCB language, these states are 
expressed as ^j(|10) — |01)) and -^(|10) + |01)), respectively. Now, the operator e~ s 
can be expressed as 



e -s _ e 9Ei=i, 2 K-i)( a *- a !) 



IT 

i=l,2 



n 



ni x t + (i - m)xl 



(B.2) 
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where Xj = ei^ - "!). Using the above, we obtain 

1 =(|10) ± |01))|m l5 m 2 ) pft = [X 1 Xl\10)±X 2 Xl\01)]\m 1 ,m 2 ) ph . 



vT 



(B.3) 



where m\ and m 2 correspond to phonon occupation numbers at site 1 and site 2 
respectively. Therefore, from equation (IB. II) we can write the density matrix element 
between singlet and triplet states in the original frame of reference as 

| «10| - <01|) ^(t) (|10> + |01» (B.4) 
= \H vh{mum 2 \ ((lOlX^I - (OllXiXl) p° T (t) (XxXjllO) + X 2 Xl|01)) \m u m 2 ) ph . 



mi ,m,2 



(B.5) 

Depending upon the presence or absence of HCB, appropriate deformation will be 
produced at each site and (XiX]|10) ± X 2 Xj|01) j \mi,m 2 ) p h represents polaronic 
states. Furthermore, in equation flB.5p . no loss in the off-diagonal matrix element on 
the left hand side implies no loss in the off-diagonal matrix element on the right hand 
side (i.e., no loss in the measured density matrix elements in the original frame of ref- 
erence) which in turn means no decoherence results. 



* Contributed equally to this work. 
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